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ABSTRACT 

This study reports pulse variation analysis results for the forth discovered 
accretion-powered millisecond pulsar XTE J1807-294 during its 2003 outburst 
observed by Rossi X-ray Timing Explorer. The pulsation is significantly de- 
tected only in the first ~90d out of ~150d observations. The pulse phase vari- 
ation is too complex to be described as an orbital motion plus a simple polyno- 
mial model. The precise orbital parameters with Porh = 40.073601(8) min and 
a^sinz = 4.823(5) It-ms were obtained after applying the trend removal to the 
daily observed 150s segments pulse phases folded with a constant spin frequency 
without Keplerian orbit included. The binary barycenter corrected pulse phases 
show smooth evolution and clear negative phase shifts coincident with the flares 
seen on the light curve and the enhancements of fractional pulse amplitude. The 
non-fiare pulse phases for the first ~60d data are well described as a fourth order 
polynomial implying that the neutron star was spun-up during the first ~60d 
with a rate z> = (1.7 ± 0.3) x 10~^^ Hz/s at the beginning of the outburst. Sig- 
nificant soft phase lags up to ~500 ^is (~10% cycle) between 2 to 20 keV were 
detected for the nonfiare pulse phases. We conclude that the anomalous phase 
shifts are unlikely due to the accretion torque but could result from the "hot 
spot" moving on the surface of neutron star. 

Subject headings: accretion, accretion disks — binaries: close — pulsars: individ- 
ual (XTE J1807-294) — stars: neutron — X-rays: binaries 



1. Introduction 



Low Mass X-ray Binaries (LMXBs) are con sidered to be the progenitor of the m il- 
lisecond pulsars detected in the radio band (see iBhattacharya fc van der Heuvell Il99ll for 
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an extensive review). It is widely believed that the weakly magnetized, slowly rotating 
neutron star is gradually spun-up through the transfer of angular momentum carried b y 
the matter from an accretion disk (lAlpar et al.lll982l : iRadhakrishnan fc Shrinivasarull982l ). 
Millisecond-time-scale variations have been observe d in many LMXB systems, including 
kilohertz Quasi-Periodic Oscillation (k Hz QPO, see Ivan der KlisI |2005| ) and burst oscilla- 



tion (see 



Strohmayer fc BildstenI |2005| for a review); however, the attempts to measure 



coherent millisecond pulsation from an LMXB, which provides a direct evolutionary link 
between radio millisecond pulsar and LMXB, were unsuccessful until the discovery of the 
first accretion-powered millisecond pulsar SAX 1808.4-3658 with a spin period Pgpin — 2.5 
ms ( jWijnands &: van der Klislll998l : IChakrabarty fc Morgan! Il998l ). in 1998. To date, coher- 



ent millisecond pu lsations have been detected in ten LMXBs with spin periods ranging from 
1.67 ms to 5.5 ms 



l isations nave been detected m ten ■LM.A..DS witn spin periods ranging iroi 
Jwiinandsll2007l : iMarkwardt et al.ll2007l : ICavriil et alTlboOTl : ICasella et al 



20071 : lAltamirano et all 120071 ) 



The fourth known accretion-powered millisecond pulsar, XTE J1807-294, was discovered 
on February 21, 2003 while the Proportional Counter Array (PCA ) onboard the Rossi X-ra y 
Timing Explorer (RXTE) monitored the Ga lactic-center region (IMarkwardt et al.l l2003al ). 
Its preliminary spin fr equency of 190.625 Hz (Markwardt et al.ll2003al ) and orbital period of 
40.0741 ±0.0005 min (IMarkwardt et al.ll2003bl ) were immediately revealed right after its dis- 



covery. Follow-up observations were made with many X-ray observatories. IMarkwardt et al. 



(l2003bh determined the source location, R. A =18^06"^59^8 pEC= - 29-24 20 (equinox 2000, 
^ 1 error) by using C/j.an6?ra observation. Campana et al. (2003), Kirsch et al.l (2004) and 
Falanga et al.l (l2005al ) reported more precise orbital par ameters and pulsation period using 
XMM-Newton observations. However, iMarkwardtl ( 20041 ) found wide swings in the apparent 
spin change rate vs. accretion rate, which were interpreted as the hot spot moving on the 
neutron star rather than from accretion torque. The wide-band (0.5 to 200 keV) spectrum 
obtained by the combination of simultaneous observations from XMM-Newton, RXTE and 
INTE GRAL is well describe d by an absorbed disk blackbody plus thermal Comptonization 
model (IFalanga et al.ll2005al ). RXTE has processed extensive monitoring of XTE J1807-294 
since February 27, 2003 with a time span of more than 150d until the source went into 
its quiescent state. The light curve exhibited an exponential-like decay with a time-scale of 
~ 19d at th e beginning of the out burst followed by a much slower decay with a time constant 
of ~ 122d (IFalanga et al.ll2005al ). Several broad soft flares lasting for hours to days were 
seen during the flux decay. The pulsed fraction increased and the pulse profile became more 
sinusoidal during the flares (jzhang et aL 2006). Twin kHz QPOs were also detected in XTE 
J1807-294 by RXTE dLinares et al.lbood l 



Owing to their fast rotation, the related parameters of the accretion-powered millisecond 
pulsars, such as orbital and spin parameters, can be precisely measured during their outbursts 
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las ting for tens to hundre d s of days. Unfortuna tely, the orbital and spin pa rameters derived 
by ICampana et~aD J2003h . iKirsch et all tooi ) and IPalanga et aP J2005ah for XTE J1807- 
294 form XMM-Newton observations were base d on the presumpt i on that the orbital period 
is fixed to the preliminary value proposed by iMarkwardt et al.l (l2003bl ) because the time 
base line is too short to constrain the orbital period. On the other hand, the ~150d's 
RXTE observations allow further binary parameter refinement, that is essential for further 
investigation. In this paper, we report the pulsation analyses of XTE J1807-294 using 
archived RXTE data (section [3]). Due to the complexity of pulse phase variation, a trend 
removal technique is employed to reveal precise orbital parameters for the system. Pulse 
phases (binary barycenter-corrected) are found to exhibit anomalous negative shifts during 
the soft fiares in all energy bands between 2-20 keV. The evolution of the non-fiare pulse 
phases shows that the neutron star was spun- up during the outburst. Significant ~ 500 /is 
soft lags for the non-fiare pulses is detected at energies from 2 to 20 keV. Analysis results, 
including possible implications for coincident soft fiares and anomalous negative phase shifts 
are discussed in section HJ 



2. RXTE Observations 

Extensive follow-up observations for XTE J1807-294 were conducted by RXTE soon 
after the source was discovered until it went into its quiescent state. We analyzed the 
available archived RXTE daXa. observed from February 27 to July 29, 2003. Th e data used to 



study the timing properties of XTE J1807-294 were collected by RXTE PCA (jjahoda et al. 



19961 ). which consists of five gas- filled Proportional Counter Units (PCUs) with the total 
collecting area ~6500 cm? and sensitive in the photon energy range 2 to 60 keV. The PCA 
data of all the observations were collected in the GoodXenon mode with a time resolution of 
~1 yus. In addition, to compare the pulsation behavior and light curve (see section [3^ . the 
data from both PCA and High Energy X-ray Timing Experiment (HEXTE) were used for 
spectral fitting to obtain the fiuxes for the all available XTE J1807-294 RXTE observations. 



3. Data Analysis and Results 

3.1. Binary Orbital Parameters 

First, all the events collected by two Event Analyzers (EAs) were combined; then, 
the event arrival times were further corrected to the barycenter of the solar system, us- 
ing JPL DE200 ephemeris and the source position determined by Chandra observation 
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( iMarkwardt et al.l l2003bl ). Only the events detected by the top layer of each PCU and 
with photon energies between 2 to 10 keV were selected in order to improve the sensitiv- 
ity. Subsequently, all the selected events were divided into 150s data segments for further 
analysis. 

The epoch folding search technique with a search resolution of 4 x 10~^ sec was applied 
to test the significance of pulsation for each 150s data segment and to obtain the preliminary 
binary orbital parameters. The test results show that the 190 Hz pulsation can be detected 
only in the first ~90d (up to TJD 12783) out of the overall ~150d's observations. In order 
to correctly estimate the orbital parameters, only those data segments with the maximum 
value from epoch folding search larger than three times of the degrees of freedom (dof=31 
for 32 bins/period, false alarm probability < 10~^) were taken into account. The trial 
period with maximum was considered the best pulsar period for each data s e gment . 



and the corresponding error was evaluated using the equation proposed by iLeahyl (119871 ) . 
The Doppler effect owing to the orbital motion can be clearly observed with the best-fitted 
mean pulsar frequency v =190.62353(3) Hz (or Ps = 5.2459421(9) ms) and binary orbital 
parameters: Porb =40.0735(2) min, sin z = 4.83(6) It-ms and =TJD12720.68245(6). 
Neither orbital eccentricity nor the spin frequency derivative can be significantly detected. 
Improved orbital parameters and the evolution of the pulsar spin frequency can be obtained 
from the pulse arrival time delay; hence the above results were applied as the initial guess 
values for the following phase analysis. 

Irrespective of whether the complex spin change could be a result of the hot spot moving 
on the surface of the neutron star or accretion torque, we can expect the intrinsic phase (after 
binary barycenter correction) to be very complicated. Therefore, we employed an alternative 
method to obtai n the precise orbital parame t ers in stead of using a method analogous to the 



one propose d bv iDeeter. Pravdq fc BoyntonI (jl98ll ) which is applied to accretion millisecond 



pulsars (e.g. iBurderi et al.ll2006l ). Because the radial velocity amplitude of the neutron star 
(K) in XTE J1807-294 is small {K/c ~ 2 x 10"^), the pulsation can be still detected in 
a 150-sec data segment if the folding frequency is close to the true spin frequency of the 
neutron star. The pulse arrival time delay due to the orbital motion can be obtained from 
the pulse phase shift. The pulse profile folded with the mean pulsar period derived from the 
orbital Doppler effect is better fitted with a two-Fourier component model (fundamental -|- 
first overtone). Typical pulse profiles are shown in Figure [H The peak of the fitted profile 
was chosen as the fiducial point for the following analysis. The error of the phase for each 
150s data segment was estimated from 1000 runs of Monte Carlo simulations. In addition to 
the time delay due to the orbital motion, significant linear phase drift (more than 1 cycle per 
day) was detected in ~3.5 hours' exposure because of slight inconsistency between the folding 
period and the neutron star spin period. To minimize the possible cycle count ambiguity. 
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we adjusted the folding period to Pjoid = 5.245942726 ms (or ujoid =190.6235070 Hz), based 
on the best pulsar ephemeris derived by the observation on TJD12700, and applied it to the 
entire data. To test if the intrinsic phase evolution is indeed complex, we fitted the observed 
pulse phase variation with the combined function of the time delay due to binary orbit plus 
evolution of intrinsic pulse phase modeled as a polynomial, that is 

C*W = <PT\To) + {uf^id -v){t- To) - h{t - T,f - hit - Tof + (1) 

However, no simple polynomial (up to 8th order) can yield acceptable fit (reduced ~ 4.2). 
We conclude that the evolution of intrinsic pulse phase is too complicated to be modeled as 
a simple polynomial. 

Since further study of the nature of accreting millisecond pulsar requires precise binary 
orbital ephemeris, we alternatively tried to remove the intrinsic phase trend for daily obser- 
vation. Although the evolution of intrinsic pulse phase was complicated for the first ~90d's 
of observations, the higher order terms of the intrinsic phase drift can be neglected for the 
daily observation because of the short exposure (^3.5 hrs). Therefore, we fitted a linear 
trend plus a single sinusoidal function due to the orbital motion with a fixed trial orbital 
period derived from the Doppler effect and then removed only the linear trend from the 
phases. The trend-removed phases were fitted with a circular orbit model to obtain updated 
orbital parameters. To further improve the orbital parameters, the new orbital period was 
used as the trial orbital period and above processes were repeated. The iterations of this 
procedure continued until improvements of all parameters were at least an order smaller than 
the corresponding errors. The final trend-removed phases are well fitted with a circular orbit 
whose best parameters are listed in Table [H The F-test shows that adding the eccentricity to 
the orbital model dose not significantly improve the fit. The F value is 1.3, which indicates 
that probability for improvement of the fit occurring by chance can be as high as ~30%. 
The pulse arrival time delay due to the binary orbital motion is shown in Figure [2l 

3.2. Pulse and Flux Correlation 

The pulse phase or spin frequency variation of the neutron star can be studies on basis 
of precise orbital parameters. For obtaining the binary barycenter-corrected pulse phases 
(hereafter, pulse phases), we folded the event time ti with a time- variable frequency due to 
the orbital Doppler effect using the parameters (Port, o-x sini, and r^/2) derived in section [371] 
as 
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where frac is the fraction of cycle counts calculated inside the braces, is folded frequency 
(see section I3.ip and Tq is the phase zero epoch of the pulsation. This operation is equivalent 
to first correcting the event times to the binary system barycenter and then folding with a 
constant frequency uq. The pulse profile for each observation ID was constructed by through 
binning the event phases. Figure [3] shows the pulse phase variation during the 2003 outburst. 
The pulse phases change like a smooth curve with several anomalous negative phase shifts 
occurring around March 9 (TJD 12707), 15 (TJD 12713), 26 (TJD 12724), 30 (TJD 12728) 
and also around May 2(TJD 12761) with time scale of several days but become somewhat 
irregular after May 11 (TJD 12770) (see top panel of Figure [3]). The first four anomalous 
phase shifts are coincid ent with the "puny " flares discovered by IZhang et al.l (120061) (see 
Figure 1 of their paper). IZhang et al.l (120061 ) also reported that the pulsed fraction increases 
with the flares. Therefore, we examined the fractional root mean square (RMS) amplitude 
for the first ~90d's observation. Figure H] shows two additional significant enhancements in 
pulsed fraction around TJD 12761 and 12769. The former one is coincident with the fifth 
anomalous negative phase shift. 

To further investigate the coincidence of anomalous phase shifts and fiares, we made a 
2-10 keV light curve for ~150d observations through spectral fitting. We used the archived 
PCA standard-2 and HEXTE data and adopted the spectral model of absorbed blackbody 
plus a power law to obtain the 2-10 keV fiuxes. This model gives reasonable fit results for all 
observations with reduced of 0.6-1.6. The bolometric correction estimated by the mean 
ratio of the integrated 2-200 keV to 2-10 keV is 2.41 ± 0.27. Figure [3] shows the 2-10 keV 
fiux of XTE J1807-294 during 2003 outburst. Four strong fiares with peak fiuxes detected 
on TJD 12707, 12713, 12724 and 12728, as well as two weak fiares with peak fiuxes detected 
on TJD12761 and 12769, can be clearly seen on the light curve. We labeled these fiares from 
1 to 6 according to their chronological order. Pulse phases exhibit negative phase shifts, and 
the pulsed fractions increase during the fiare states. From Figure [3] and Figure HI the fiares 
are likely to appear in pairs with peak separations of 5.9d for fiares 1 and 2, 4.8d for fiares 
3 and 4 and 7.8d for fiares 5 and 6. 

These negative phase shifts are more likely due to the motion of the accretion footprint 
rather than the spin-down torque acting on the neutron star, which will be discussed in 
section 14.21 Further analysis results for the pulse behaviors during the fiare states are 
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presented in section 13.31 On the other hand, the smooth pulse phase evolution for the non- 
flare state during the flrst ~60d of the outburst is possible owing to the accretion torque. We 
deflned the epochs for the flare states: TJD 12703 to 12716 for flares 1 and 2, TJD 12720 to 
12730 for flares 3 and 4 and TJD 12759 to 12778 for flares 5 and 6. To investigate the smooth 
pulse phase evolution during the non-flare states, we removed the phases during the flare 
states and also discarded the ones after TJD 12778. The non-flare phases were fltted with a 
polynomial, starting from a quadratic function and then subsequently adding higher order 
terms. The F-test was applied to verifled if the additional higher order term signiflcantly 
improves the fltting (with confldence level > 90%). We found that the phase evolution for 
the non-flare states is better described by a 4th order polynomial (see Figure [5]). On the 
other hand, imperfect solar system baryc enter correction due to source position error results 
in correction error for event arrival time (ILyne &: Graham-Smitlj|l998f). which further leads 
to systematic errors in spin parameters. From Eq. 5.3 in lLyne &: Graham-Smithl (Il998l ). the 
systematic errors of the pulse phase and spin parameters, including phase zero epoch (Tq), 
spin frequency (u) and its derivatives {d^i'/dt'^) can be expressed as 



IF 



:^sin[fi^(t-r,)-A-/3] 



STn = 



V c 

d{5(t)) _ 
~dt 
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where Re is the orbital radius of the Earth, fi^; is Earth's orbital angular velocity, and is 
the time of vernal point. The best fltted spin parameters and the corresponding statistic and 
systematic errors are listed in Table [2l The frequency derivative i>{t) derived from the best 
fltted spin parameters shows that the neutron star was spun-up during the flrst ~60d of the 
2003 outburst (see Figure [5]). Further, the non-flare light curve can be well described as the 
combination of a fast and a slow exponential decay component with decay time constants of 
Tf = 11.4 ± 0.2 d and = 95 ± 2 d respectively (see Figure [3]). 
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3.3. Energy-dependent Pulse Behaviors 



3.3.1. Soft Lags for Non- flare Pulses 



Cui. Morgan fc TitarchukI (119981 ) discovered that a high-energy pulse leads a low-energy 
pulse by about 200 /is (referred to as "soft lags") from 2 to 10 keV but saturates at a higher 
energy for the accreting millisecond pulsar SAX J1808. 8-3658. T his energy-dependent pulse 



arriva l time has been well explained by t he scenario proposed by iGierlihski. Done fc Barret 



(I2OO2I ) and lPoutanen &: Gierlihski The emissions mainly consist of a soft blackbody 

component from a heated spot on the neutron star surface and a hard component from 
the Comptonization of the blackbody photons (seed photons) in the plasma heated by an 
accretion shock slab with Thomson optical depth of 0.3 to 1. The hard Comptonization com- 
ponent leads the soft blackbody component as a result of the different angular distribution 
of the soft and hard components and the Doppler boosting effect. The low-energy pulses lag 
behind the high-energy pulses because the latter contain relative less soft blackbody compo- 
nent. For the photons energies greater than ~10 keV, only the Comptonization component 
is left so the hard leads saturate. The emission spectra of many accretion-powered mil- 
l isecond pulsars are w ell described by this two-component model, including XTE J1807-294 
(IFalanga et al.ll2005al ). Subsequent to the discovery of soft lags in SAX J1808. 8-3658, similar 
phenomenon was also seen in some of accreting millisec ond pulsars, such as X TE J1715-305 
Jcierhhski fc Poutanen|[2005h and HETE J1900. 1-2455 JCallowav et~aDl2007h . However, for 
IGR J00291-I-5934, the soft l ags increase up to 7 keV and then decrease at higher energy 
bands instead of saturation (IFalanga et al.l l2005bl ) . The decre ase at higher energy bands 



is likely due to the t ime lags caused by Compton upscattering (ISunyaev fc TitarchukI Il980 
Vaughan et al.l 119971 ). 



To investigate the energy-dependent pulse behaviors for XTE J1807-294, we divided 
2-20 keV photons into 8 energy bands for the first ~60d observations. Since we assumed 
hot spot to be fixed during the non-fiare state, we folded the non-fiare event times of each 
energy band with a ephemeris (hereafter, the best ephemeris) from the variable observed 
pulse frequencies due to the orbital Doppler effect and neutron star spin parameters from 
the 4th order polynomial fitting (see Table [2]). Thus, pulse phases of the non-fiare state are 
expected to be a constant with time. The average pulse profile for the non-fiare state of each 
energy band was obtained by binning the corresponding event phases. The energy-dependent 
pulse arrival time delays, measured relative to the pulsation of the softest band (2-3.6 keV), 
were obtained through the cross correlation of the pulse profiles. The high energy pulses lead 
soft pulses up to ~ 500 ^is (~ 0.1 cycle) from 2 to 10 keV and the leads saturate for 10 to 20 
keV pulse (Figure [HD, similar to that seen in SAX J1808. 8-3658 (ICui. Morgan &: Titarchuk 
19981 ). On the other hand, the fractional RMS pulse amplitude slightly increases from 4.5% 
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to 5% in the 2-4 keV range, remains constant to 7.5 keV, and largely decreases to only 2% 
at 20 keV (see Figure [7]). 



3.3.2. Phase Shifts during Flare States 

Significant phase shift, and its correlation to the flux and fractional pulse amplitude, 
has never observed for other accreting millisecond pulsars. To further understand its origin, 
we studied the energy-dependent phase shifts during the flare state. The pulse phase defined 
here is identical to the one in section 13.3.11 that is, folded with the best ephemeris. The 
phase shift for each energy band, defined as the phase difference between the flare state 
and the corresponding mean phase of the non-fiare state, was obtained through the cross 
correlation method. However, since this ephemeris is applicable for the first ~60d data, only 
the phase shifts of the first four flares can be well-defined. Figure [8] shows that the phase 
shifts of different energy bands when the fluxes of flares attained their maximum values for 
the flares 1 to 4. Significant phase shifts can be observed in all energy bands for all the four 
flares. In general, the phase shifts decrease for higher energy bands from 2 to 5 keV, almost 
a constant between 5-11 keV and slightly increase for the hardest energy band. 



Discussion 



4.1. Orbital Parameters 



We have derived the precise binary orbital parameters for th e X-ray binary XT E J1807 



294. The orbita 



para n ieters have b e en pr evious ly reported by iMarkwardt et al.l (l2003al ) 



Campana et al.l (120031 ) , iKirsch et al.l (120041 ) , and iFalanga et al.l (l2005al ) , whose values are 



listed in Table [3] for comparison. All the parameters are consistent with each other except 
the epoch of the mean longitudes for the binary orbit (i.e. phase zero epochs). The epoch of 
mean l ongitude is chosen as when mean longitude=90° (Tt,/2) in this paper. iFalanga et al. 



(l2005al ) claimed that they used same reference point, but it deviates from ours by about 
~0.5 phase. We examined their arguments and found that t hey probably u s ed the positive- 
crossing time in time vs. pulsar period plot (see Figure 6 in IFalanga et al.l (l2005al )). which 
is equivalent to the mean longitude=270°, as the phase zero epoch. If their reference time 
is indeed T37r/2 instead of , the value proposed by them and the one in t his paper are 



consistent. On the other hand, the epoch of mean longitude considered by iKirsch et al. 



(120041 ) is equivalent to the mean longitude=0°, which is close but still signif icantly different 



from what we found (more than 8a). Their value, and the one proposed by ICampana et al. 
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( 20031) ■ is probably due to a short observation time span and, as mentioned in iFalanga et al. 
( l2005al ). the complex method used by these authors. 



4.2. Phase Shift during Flare State 

With the precise orbital parameters, we have found significant pulse phase shifts on the 
time scale of one to several days coincident with the fiares as well as pulsed fraction enhance- 
ments. If the hot spot is assumed to be fixed on the neutron star surface and the pulse phase 
variation truly refiects the history of the change in the spin frequency of the neutron star, fit- 
ting a quadratic curve for the pulse phases around the fiare states implied a spin-down rates 
10~^^ Hz/s. It is much larger th an the theoretical expectation value (~ 10~^^ Hz/s, see 
Rappaport. Fregeau fc Spruitll2004j ) or t he reported values fo r accretion-powered millisecond 
pulsars (~ 10~^^ — 10~^^ Hz/s, see e.g. iBurderi et al.l 120071 ). In addition, if the fiares are 
caused by an enhanced accretion, the ad ditional accretion torque should further spin-up the 



neutr on star rather than spin it down (iGosh fc LambI Il979l : iRappaport. Fregeau fc Spruit 



20041 ). Thus, the fiare- coincident negative phase shifts are highly unlikely due to the accre- 
tion torque. 

Another possible explanation for the anomalous pulse phase shifts is the hot spot moving 
on the surface of neutron star. Assume that the neutron star is rotating eastward as the 
Earth and that the hot spot is fixed at a position in the northern hemisphere during the 
non-fiare state. As the accretion rate increases in a fiare state, the hot spot moves eastward, 
in the same direction as that of the neutron star rotation, so that the pulse peak arrives to 
an observer earlier than the one expected for the non-fiare state and results negative phase 
shift. At the time when the fiare reaches its maximum, this movement can be as large as 
about 70° (~0.2 cycle) in longitude. When the fiare's fiux decreases, the hot spot moves 
westward and eventually returns to its original location at the end of the fiare. 

The correlation between phase shift and fiux has been suggested for other accreting 
millisecond pulsars, although the no fiare similar to tha t observed in XTE J 1807-294 has 
ever been detected in other accreting millisecond pulsars. iPapitto et al.l (120071 ) reported the 
discovery of high anti-correlation between X-ray fiux and the pulse phase jitter (~ ±0.05) 
around the mean quadratic phase trend in the accreting millisecond pulsar XTE J1814-338. 
Such a correlation is very similar to what we found in XTE J1807-294 during the fiare states 
although the correlati on between th e pulse d fraction and the fiux is still unknown for XTE 
J1814-338. Recently, iBurderi et al.l (120061 ) discovered an anomalous large phase jump for 
SAX 1808.4-3658 during its 2002 outburst, which is also likely as a result of a large accretion 
rate change (see Figure 1 in IBurderi et al.ll2006l ). 
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Zhang et al.l (120061 ) suggested that the flares seen in XTE J1807-294 could result from 
accreting blobs of high mass density, lasting for hours to days, at the inner edge of the 
accretion disk. The positive correlation of kHz QPO frequency and pulsed fraction implies 
that the s urface emission pro perties of neutron star are very sensitive to the inhomogeneous 
disk flow (IZhang et al.l l2006l ) . We found that not only the pulsed emissions increase but 
also the position of the hot spot moving on neutron star surface. Since the longitudinal 
movement of the hot spot on the neutron star results in the phase shift, it is also possible 
that the hot spot moves in latitude, which could be a part of the reason for the pulse 
amplitude enhancements during the flare states. 

The flux and ph ase shift coincidence ha s also been seen in accreting pulsars with much 
lower spin frequency. IStrickman et al.l (119961 ) reported that during the bursting X-ray pulsar 
GRO J1744-28 {u = 2.1Hz) type II bursts, likely a consequence of accretion instability, the 
pulse arrival time lags by ~ 90 ms relative to the non-burst state. Pulsed fraction is also 



strongly enhanced during the bursts (IKouveliotou et al.lll996l ). iMillerl (119961 ) proposed a 
scenario to describe the pulse phase shifts. Because of the misalignment of the neutron star 
spin axis and magnetic dipole moment, the accretion flow is asymmetric and its footprint 
forms an arc surrounding the magnetic pole. The enhanced accretion causes the footprint 
to move around magnetic pole as a consequence of the deformation of the field lines. The 
azimuth movement of the footprint with respect to the spin axis results in the observed phase 
shift. Except that the phase shift is negative for XTE J1807-294 but positive (pulse arrival 
delay) for GRO J 1744-28 and the different time scales for flare (several days) and type II 
burst (several seconds), the coincidence of phase shift, pulse amplitude enhancement and 
burst in GR O J1744-28 is very similar to what we report for XTE J1807-294. The scenario 
suggested by iMillerl (119961 ) can probably give us a clue to explain the observed behaviors of 
XTE J1807-294. 

Form the discovery of the coincidence between the flare and phase shift in XTE J1807- 
294, as wel l as the correlatiori between the X-ray flux and the pulse phase jitters in XTE 
J1814-338 (jPapitto et al.ll2007l ). hot spot posi t ion on accreting millisecond pulsar is sensitive 
to accretion flow, as suggested by lZhang et al.l (120061 ). Although we assume that the accretion 
hot spot is still fixed on the surface neutron star during the non-flare state, we cannot exclude 
the possibility that the hot spot marginally drifts as the X-ray flux changes. Therefore, it 
is difficult to determine whether the pulse phase variation is a consequence of the spin 
frequency change or hot spot drift on the neutron star surface. If the correlation between 
the hot spot drift and accretion flow is applicable to all accreting millisecond pulsars, the 
measurements of the spin parameters (e.g. z> ) from which accretion torque is estimated, are 
not reliable, including the ones listed in Table [2] for XTE J1807-294. This may also explain 
the large disparity between the measured spin frequency derivative and the theoretical value 



calculated form the accretion torque reported by lMorgan et al.l (120031 ) for SAX J1808. 4-3658 
in its 2002 outburst. 



5. Conclusion 



1. Due to complicated intrinsic phase variations, to obtain the precise orbital parameter, 
we folded our data sets with a constant frequency and used the trend removal technique 
to remove the local linear pulse phase trend from the observed pulse phase variation for 
daily observations. We obtained more precise orbital parameters consistent with the ones 
previously reported (except for T^/2). 



2. The binary barycenter-corrected pulse phases change smoothly except for several anoma- 
lous negative phase shifts coincident with flux flares and pulse amplitude enhancements. The 
pulse fractional amplitude variation indicates that there were six flares during the first ~90d 
of the XTE J1807-294 in its 2003 outburst. The phase changes of non-fiare pulse imply that 
the neutron star was spun-up during the first 60d of the outburst under the assumption that 
the hot spot was fixed on the surface of the neutron star during the non-fiare states. 



3. The analysis results of the energy- dependent pulse arrival time indicates that hard 
pulses lead soft pulses up to ~ 500 fis from 2 to 10 k eV, and the lead saturate s frorn . 
10 t o 20 keV, which agrees with t he scenario proposed by iGierlihski. Done &: Barretl ((2002) 
and iPoutanen fc Gierlihskil (120031 ) . 



4. Our discovery of the coinci dence of flare and phase shift, as well as its correlation with 
pulse amplitude and KHz QPO Jzhang et al.lbooeh in XTE J1807-294, implies that probably 
due to weak magnetic confinement on the accreting material, the hot spot position on the 
accreting millisecond pulsars is sensitive to the accretion flow. 
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Fig. 1.— The pulse profiles (2-10 keV) of XTE J1807-294 measured on TJD12697 (upper 
left), TJD12712 (upper right), TJD12720 (lower left) and TJD12775 (lower right). The solid 
line on each plot is the pulse profile fitted by the two-Fourier component model. 
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Fig. 2. — Pulse arrival time delay (after the phase trend removal) due to the orbital motion 
(top) and the fit residuals (bottom) . 



- 18 - 




50 100 150 

TJD- 12700 



Fig. 3. — Comparison between evolution of 2-10 keV pulse phase folded with a constant 
frequency plus Keplerian orbit (i.e. Eq[2]) (top panel) and 2-10 keV flux (bottom panel). A 
total six flares were detected during the flux declination. The arrows in both panels indicate 
the flare peak time. Clear flare and phase shift coincidence can be seen. The solid line on the 
bottom plot represents the best fitted two-component exponential model for the non-flare 
state fluxes. 



- 19 - 




-20 -10 10 20 30 40 50 60 70 80 90 100 

TJD-1 2700 



Fig. 4. — RMS pulsed fraction vs. time. Similar to Figure [3l the arrows indicate the times 
for the six flare peaks, showing the pulse amplitude and flare coincidence. 
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Fig. 5.— Top: Phase evolution of non-flare state from TJD12697 to 12757. The sohd hne 
is the best-flt 4th order polynomial. Bottom: The spin frequency derivative {dotv) derived 
from the best-fit parameters from phase evolution, assuming that the hot spot is fixed on 
the surface of neutron star. 




Fig. 6. — Energy- dependent pulse arrival time delay relative to the softest energy band (2-3.7 
keV) for the non-flare pulse. Negative value indicates that the pulse leads the softest one. 
The leading saturates at ~10 keV about 500 /xs (~0.1 cycle). 
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Fig. 7. — Energy- dependent mean pulsed fraction for the non-flare state. 
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Fig. 8. — Phase shifts of different energy bands for the flare states measured on the first four 
flare peak times. The flare index is labeled on the up-left corner of each plot. Signiflcant 
phase shifts can be seen in higher energy bands, although the lowest energy band has the 
largest phase shift 



Table 1: Orbital Parameters for XTE J1807-294 from RXTE Observation 



Parameter Value 

Orbital period, Porb (min) 40.073601(8) 

Projected semimajor axis, sin i (It-ms) 4.823(5) 

Epoch of 90° mean longitude, (TJD/TDB) 12,720.682574(6) 

Mass function (M©) 1.497(4) x 10"^ 

Eccentricity, e <0.006(2f7) 

xVDegree of freedom 2171/1880 



- 25 - 



Table 2: Neutron Star Spin Parameters for XTE J1807-294 



Parameter 


Value 


Statistic 


Systematic 






Error 




Error ^ 


Phase zero epoch, Tq (TJD/TDB) 


12700.00000000(2) 


1.0 X 10" 


^10 


1.7 X 10-^ 


Frequency at Tq, z/q (Hz) 


190.62350686(8) 


1.4 X 10- 


-8 


7.4 X 10-^ 


dv/dt at Tq (Hz/s) 


1.7(3) X 10^13 


3.0 X 10- 


-14 


1.1 X 10-1^ 


d'^v/de at To (Hz/s^) 


-1.1(1) X 10-19 


1.0 X 10' 


-20 


2.9 X 10^21 


d^v/dt^ (Hz/s^) 


4.4(4) X 10"26 


4.0 X 10" 


-27 


4.3 X 10-2^ 


/Degree of freedom 


22.76/26 









"The parameter errors are evaluated from quadratic sum of statistical and systematic errors 

, // 

The systematic errors are estimated from the ~ 1 source position error (see 13. 2p 
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Table 3: Measured Orbital Parameters for XTE J1807-294 



Mission 


p 

orb 


a sinz 


Epoch of mean 


Phase for epoch 


Reference 




(min) 


(It-ms) 


longitude (TJD) 


of mean longitude^ 




RXTE 


40.0741(5) 








1 


XMM-Newton 


40.0741(fixed) 


4 7+1-2 

^- ' -0.04 


12720.68644(13) 


0.140(5) 


2 


XMM-Netwon 


40.0741(fixed) 


4.8(1) 


12720.67415(16) 


0.699(6) 


3 


XMM-Netwon 


40.0741(fixed) 


4.75(39) 


12720.668(2) 


0.48(7) 


4 


RXTE 


40.073601(8) 


4.823(5) 


12720.682574(6) 




5 



"Folded with the orbital ephcmcris proposed in this paper. See l4.1l for the discussions about the inconsistencies 
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